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Abstract. We investigate transport in a three-terminal graphene quantum dot. All 
nine elements of the conductance matrix have been independently measured. In the 
Coulomb blockade regime accurate measurements of individual conductance resonances 
reveal slightly different resonance energies depending on which pair of leads is used for 
probing. Rapid changes in the tunneling coupling between the leads and the dot due 
to localized states in the constrictions has been excluded by tuning the difference 
in resonance energies using in-plane gates which couple preferentially to individual 
constrictions. The interpretation of the different resonance energies is then based on 
the presence of a number of levels in the dot with an energy spacing of the order of 
the measurement temperature. In this mult i- level transport regime the three-terminal 
device offers the opportunity to sense if the individual levels couple with different 
strengths to the different leads. This in turn gives qualitative insight into the spatial 
profile of the corresponding quantum dot wave functions. 
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1. Introduction 

Craphene nanostructures are believed to have potential applications in both 
conventional electronics and solidstate quantum information processing. In particular, 
graphene quantum dots might be promising for spin qubits due to their predicted long 
spin life times. [IJ 

As a consequence of the gapless band structure charge carriers cannot be 
electrostatically confined in graphene. However, by cutting graphene into narrow 
ribbons a so-called transport gap is opened where the current is suppressed around 
the charge neutrality point. [2[ [3l |4[ [5l El [Tj By using short and narrow constrictions as 
tunnel barriers more complicated nanodevices like quantum dots have been successfully 
created. This has led to a number of experiments where, for example, excited states have 
been observed in single ^ and double quantum dots [9l [THl [H] , spin states have been 
investigated [12j and the electron-hole crossover has been studied [13]. In addition, the 
modulation of transport through graphene quantum dots due to localized states in the 
constrictions have been investigated in several studies. [TH [151 US] Still, there are open 
questions concerning the detailed influence of constriction localizations on transport on 
small energy scales. 

The current through a two-terminal quantum dot does not give access to the 
individual coupling strengths between the dot and each lead. However, if a dot in 
the single-level tunneling regime of the Coulomb blockade is connected to three or more 
leads the individual tunnel coupling constants between each lead and the quantum dot 
can be determined from measurements of the conductance matrix of the system. [IT] 

Following the approach of Ref . [17] , we here investigate transport in a three-terminal 
graphene quantum dot in the multi-level regime. The three terminals offer the possibility 
for fast and convenient probing of the conductances of each lead, thereby providing 
further insight into how localized states in the constrictions affect transport through 



CONTENTS 



3 



V, 



bias 3 



9 



3 



1 




1 



V, 



bias 1 




bias 2 



1 



Figure 1. (a) Scanning force micrograph image of the measured quantum dot with a 
sketch of the measurement setup. Leads 1, 2 and 3 are highhghted in red, orange and 
black, respectively. Three plunger gates are used to tune the dot in addition to the 
global back gate. To each lead a bias voltage can be applied and the current flowing 
can be measured. 



the dot. In addition, being in the multi-level regime gives us the unique chance to 
experimentally observe how different leads couple with different strengths to different 
dot states. 

2. Sample and experimental methods 

Single layer graphene flakes were exfoliated from natural graphite, deposited onto a 
highly doped Sihcon substrate covered by 285 nm thermal silicon dioxide, and identified 
using Raman spectroscopy [181 HH] ^^id light microscopy. In a first electron beam 
lithography (EBL) step, followed by metal deposition of 5nm titanium, 45 nm gold and 
lift-off, the Ohmic contacts were added to the flake. The structure is then patterned 
by a second EBL step followed by reactive ion etching with argon and oxygen (for a 
detailed description of similar fabrication see Ref. [20j). 

A scanning force micrograph (SFM) image of the measured quantum dot is depicted 
in Fig.[TJ The quantum dot is connected to three leads, labeled 1, 2 and 3, through 
narrow constrictions. From the SFM image the diameter of the dot is determined to 
be 110 nm and the width of the constrictions is found to be 40 nm. In addition to the 
global silicon back gate (BG) three in-plane plunger gates, PGi, PG2 and PG3, are used 
to tune the dot and the constrictions. The remaining three in-plane gates influence the 
transport through the dot only weakly and are therefore not used. 

In Fig.[T] we additionally sketch the measurement setup. In all measurements 
presented in this study a DC bias voltage is applied to one of the three leads while 
the other two leads are grounded. The currents through the three leads are measured 
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Figure 2. (a) Current through the three leads as a function of back gate voltage for 
a large backgate voltage range. A 1 mV bias voltage is applied to lead 1, while leads 
2 and 3 were grounded. A transport gap of ~ 12 V in back gate can be seen, (b)-(d) 
I3, I2 and Ii as a function of Vbg and Vpgi for large gate voltage ranges. The broad 
diagonal lines are due to resonances in the constrictions, (e) A zoom of (d) where the 
narrow diagonal lines corresponds to Coulomb peaks. 



simultaneously using current-voltage converters. All measurements are performed at 
1.7 K unless stated otherwise. 

3. Results and discussion 

3.1. Device characterization 

Fig.[2](a) shows the currents through the dot for a large range in Vbg- Curves labeled 
/i, I2 and /3 correspond to the current measured in lead 1, 2 and 3 respectively. For this 
measurement ImV bias voltage is applied to lead 1, while leads 2 and 3 are grounded. 
We use the convention that negative currents flow from the leads into the quantum dot, 
while positive currents flow from the quantum dot into the leads. 

Around the charge neutrality point a region of 12 V of suppressed current, 
corresponding to the transport gap [4j, can be seen. Within this region Coulomb 
blockade is observed. From Coulomb diamond measurements we determine the charging 
energy of the quantum dot to be 8-15 meV (not shown). 

It can also be seen that constriction 2 is generally more closed than 1 and 3. For 
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Table 1. Relative lever arms apc/ci^BG for the in-plane gates with respect to the dot 
and the three constrictions. 

^.Q^ /a^.Q^ ^Gonstr.l /^Constr.l ^ Constr.2 Gonstr.2 ^ Constr.3 /^Gonstr.3 

^PG/^BG ^PG /^BG ^PG / ^BG ^PG / ^BG 

PGi 0.59 1.15 0.68 0.25 

PG2 0.50 0.13 0.88 0.65 

PG3 0.58 0.65 0.13 1.15 



high charge carrier densities I2 is less than 10% of Ii (the total current). This asymmetry 
is also present in the regime of Coulomb blockade where the current flowing through 
constriction 2 is often too small to be measured and the quantum dot is effectively a 
two-terminal dot. Still, it is possible to find regimes where the current contributions 
from the three leads are comparable and in the following sections we will focus on one 
of these regimes. 

In order to characterize the device further we measure the current through the 
three leads as a function of V^q and each of the three plunger gates on both a small 
and a large voltage scale, and determine the plunger gate lever arms relative to the back 
gate lever arms q^pg/<^bg with respect to the dot and each of the three constrictions. 
As an example Fig.[2](b)-(d) shows /i, I2 and respectively, as a function of V^q 
and VpGs • A voltage of 1 mV is applied to lead 1 and Ii , I2 and /a are measured 
simultaneously. The gates are swept over a large voltage range and the broad diagonal 
lines that are visible are attributed to resonances in the constrictions. [14j In Fig.[2](d) 
three different slopes (marked by orange lines) can be identified, while only two different 
slopes can be found in Fig.[2](b) and (c). Since lead 1 is biased all charge carriers 
flowing through the dot have to tunnel through constriction 1. As a result, resonances 
originating from localized states in constriction 1 are seen in all three currents. On 
the other hand, a resonant state in lead 2 or 3 will enhance the current only in this 
speciflc lead, with the consequence that only the current through this specific lead 
(and the biased lead) is enhanced. From this measurement we can therefore assign 
one slope to states in each constriction and subsequently determine o^pg2/<^bg for all 
three constrictions. Complementary measurements were done for the two other plunger 
gates and the complete set of q^pg / Q^bg is summarized in Table[l] These lever arms are 
consistent with the geometry of the sample (Fig.[T]). 

Fig.[2](e) shows a high resolution measurement corresponding to a zoom of Fig.[2](d) 
(see the black square in Fig[2](d)). The narrow diagonal lines correspond to single 
Coulomb resonances in the quantum dot and from the slope we determine o^pg/<^bg 
for the dot. From corresponding measurements varying the other two in-plane gates the 
relative dot lever arms of all three plunger gates are extracted. These lever arms can 
also be found in Table[Tl 

From Table[T] it can be seen that the relative lever arms with respect to the dot 
^pg/^bg ^^"^y similar for all three plunger gates. However, the lever arms with 
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Figure 3. (a)-(c) Measurement of the complete conductance matrix for the system. 
The conductance in lead 1, 2 and 3 is plotted in red, orange and black respectively. 
The conductance in the biased lead is always plotted negative, (d) Corresponding 
individual conductances. 



respect to the different constrictions vary significantly. In particular the lever arm 
of each plunger gate with respect to the constriction on the opposite side of the dot is 
much weaker than all other lever arms. Hence, the plunger gate dependence of transport 
through the dot can be used to identify if changes are due to alterations of the dot wave 
function or the constriction resonances. 



3.2. Determination of individual conductances from the conductance matrix 
The conductance matrix G of a three-terminal system is given by 
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In Fig.|3[a)-(c) the nine elements of the conductance matrix are shown, measured by 
applying a 100 uV bias to lead 1, 2, and 3 respectively. There are two sum-rules 
that should be obeyed by the conductance matrix. First, due to current conservation 
Tif^^Gij = for all j. Second, if the same voltage is applied to all leads no current 
should flow, Ti^-^-^Gij = for all i. In addition, at zero magnetic field G should be 
symmetric, Gij = Gji. As a result there are only three independent conductance matrix 
elements, from which the complete matrix can be deduced. For the measurement shown 
in Fig.[3](a)-(c) the first sum rule is obeyed with a relative error less than 1 % of the 
highest current level, while the second sum rule is obeyed with a relative error less than 
10 % of the highest current level. In order to obtain such a small error the measurements 
are done very carefully. To minimize the influence of voltage offsets in the measurement 
setup, measurements for positive and negative bias were averaged. In addition, to avoid 
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errors due to instabilities of the sample all nine elements of the conductance matrix for 
both positive and negative bias are measured before each BG step. In other words, all 
conductance matrix elements are measured within a single back gate sweep. 

For single-level transport in the weak coupling regime the individual tunnel 
couplings, r, between the dot and each lead can be determined from the conductance 
matrix. [21j In this transport regime the width of the Coulomb peaks is peak- 
independent for a given temperature. Looking at the Coulomb peaks in Fig.[3](a)- 
(c) it can be seen that in our measurements the width of the peaks varies (especially 
pronounced for the forth peak). This is a sign of multi-level transport. 

In the multilevel regime the individual tunnel coupling strengths from the dot to 
each lead cannot be extracted directly from the conductance matrix. However, if we 
consider the quantum dot as a classical star-shaped conductance network, we can extract 
the three individual conductances Gk connecting lead k to the dot from the relation 



G 



Gi + G2 + G3 



^ Gi{G2 + G3) —G1G2 — G1G3 ^ 

—G2G1 G2{Gi + G3) — G2G3 
^ —GsGi — G3G2 ^3(^1 + G2) J 



(2) 



In Fig.[3|(d), Gi, G2 and Gs obtained from the nine conductance matrix elements 
shown in Fig.|3](a)-(c), can be seen. The individual conductances fluctuate largely from 
peak to peak demonstrating that the coupling strengths between the leads and the dot 
vary signiflcantly from peak to peak. 



3.3. Temperature dependence 

In the previous section signs of multi-level transport were seen. To further support this, 
we present the temperature dependence of the Coulomb peaks shown in Fig. [3) This is 
depicted in Fig. [4] where the three independent conductance matrix elements G21, Gsi 
and Gs2 are plotted as a function of Vbg for seven diflFerent temperatures between 1.7 K 
and 10 K. In general it can be seen that all Coulomb peaks broaden with increasing 
temperature. For the flrst peak (from the left) and the third peak, the peak maxima 
increase for increasing temperature, which is a signature of multilevel transport [22j. 
However, for the second peak the peak maximum increases with increasing temperature 
for Gsi and G32 while it decreases for G2i- A similar behaviour is seen for the fourth 
peak where the peak maximum increases with increasing temperature for G21 and G31 
while it decreases for G32. It should be noted that even though two peaks are seen to 
decrease in height, they do not show the 1/T-dependence as expected for true single-level 
transport. 

It is known that in the multi-level regime the temperature dependence of Coulomb 
peaks can vary from peak to peak due to variations in the couplings between the leads 
and the diflFerent dot states. [22] However, the measurement of diflFerent temperature 
dependences of conductances measured in diflFerent leads for the same Coulomb peaks is 
unique to a three or more terminal system and has to our knowledge not been measured 



CONTENTS 



8 




Figure 4. Temperature dependence of the four studied Coulomb peaks. The three 
independent matrix elements G21, G31 and G32 are plotted. For the first and the third 
peak (from the left) the peak maxima increase with increasing temperature for all 
three conductance matrix elements. However, for the second peak the peak maxima 
increase with increasing temperature for G31 and G32 while it decreases for G21. A 
similar behaviour with different temperature dependence for the different conductance 
matrix elements are seen for the fourth peak. 



before. Following the arguments of Ref. [22j, our results suggest that the different leads 
couple with different strengths to the different dot states involved in transport. This 
is further supported by detailed measurements of single Coulomb peaks as discussed 
below. 

3.4' Shift between Coulomb resonance positions due to coupling of different leads to 
different dot states 

When studying single Coulomb peaks in detail we frequently observe that peaks 
corresponding to conductances in different leads have their maxima at slightly different 
positions in gate-voltage. An example of this is shown in Fig.[5](a) where G31 and G21 
are plotted for the fourth peak in Fig.[3](a) (see red star). A 100 //V bias voltage is 
applied to lead 1, a voltage of 100 mV is applied to gate 1 while all other gates are 
grounded. The conductance G31 (black dotted curve) has its maximum at 4.217 V while 
G21 (orange solid curve) has its maximum at 4.215 V. For multi-level transport the full 
width at half maximum (FWHM) of a Coulomb peak is 4.4 /cbT, where is the 
Boltzmann constant and T is the temperature [2T]. From the FWHM of the measured 
peaks we then estimate the shift between the maxima of G31 and G21 to be 0.9 /cbT. 

It should be noted that the measurements shown in Fig.[5](a) are not averaged 
between conductances measured for positive and negative bias (unlike in Figjs]). The 
gating effect of the lead where the bias is applied shifts the Coulomb resonances in 
energy at finite bias voltages. For positive bias voltages the resonances are shifted 
to more positive back gate voltages, while for negative bias voltages they are shifted to 
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Figure 5. (a) Measurement of G21 and G31 as a function of Vbg foi" a single 
Coulomb peak. The maxima of the two peaks are shifted by 2 mV in back gate voltage 
corresponding to 0.9 k^T when assuming that the FWHM of the peaks are ^ 4.4 k^T. 
(b) Corresponding calculation with a two-level model where the three leads couple 
differently to the two dot levels. 

more negative back gate voltages. Averaging would therefore result in a broad resonance 
with a maximum positioned between the maxima of the original resonances. By only 
considering the conductances measured with the bias voltage applied to the same lead 
(for positive or negative bias voltage), the gating effect of the source causes the same 
shift for all resonances. Hence, the shift of 2mV between the maxima of G21 and G31 
found above is not due to the gating effect of the source. 

The shift between the maxima of the conductances of G21 and G31 can be 
understood in terms of multi-level transport where the leads couple with different 
strengths to the different dot states. In order to illustrate this effect qualitatively, we 
calculate G31 and G21 for the simplest possible multi-level system, a two level system. 
We use the rate equation approach introduced by Beenakker in Ref. [2T] extended to a 
three-terminal dot with two levels contributing to the current. 

In this model it can be shown that the current in each lead is the sum of two 
contributions, one current via the first level and one current via the second level. The 
shift between currents in two different leads is determined by three parameters, the 
single particle level spacing A and two parameters determining how the current in each 
lead is distributed between the two dot levels. The measured shift can be qualitatively 
reproduced by the model for a large range of parameters. An example of a calculation 
showing good agreement with the experiment in Fig.[5](a) is depicted in Fig.[5](b). In 
order to put more constraints on the values of the parameters we also tried to reproduce 
the temperature dependence of the peaks. Unfortunately, with a two level model it is not 
possible to quantitatively reproduce the observed shift and the observed temperature 
dependence at the same time. Thus, we here most likely have more than two levels 
involved in transport. Still, we would like to emphasize that the simple two-level model 
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Figure 6. (a) Shift relative to the average FWHM of the peaks as a function of 
FpGi- The red Une is a hnear fit to the data points, (b) Corresponding evolution of 
the individual conductances as a function of Fpgi- G2 hardly changes, Gi changes 
only slightly, while G3 changes significantly. For the evolution of the individual 
conductances as a function of Fpgs (c) and l^Og (d) Gi and G2 again change only 
slightly while G3 changes significantly. 



does qualitatively reproduce the shift, supporting that it is indeed due to the different 
coupling of different leads to different dot states. This also agrees with the interpretation 
of the temperature dependence of the Coulomb peaks as discussed above. 

3.5. Evolution of Coulomb resonance shift with in-plane gate voltages 

From the above discussion we argue that the observed shift is due to different leads 
coupling with different strengths to different quantum dot states. However, it is still an 
open question if small changes in the coupling strengths between the leads and the dot 
are dominated by changes of the localized states in the constrictions or changes in the 
dot wave function. 

Fig.[6](a) shows the shift from Fig.[5](a) relative to the average FWHM of the peaks 
as a function of l/pGi- It can be seen that already for small plunger gate voltages 
the shift can be tuned significantly and systematically by an in-plane gate. Fig.[6](b) 
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shows the evolution of the individual conductances Gk of the three leads with Vy>g^. 
Fig.|6](c) and (d) show the corresponding evolution of Gi, G2 and G3 for l/pcs and Vy>g^ 
respectively. If the changes in the coupling between the dot states and the lead states 
would be dominated by changes in the localizations in the constrictions we would expect 
a correlation between the evolution of the individual conductances and the relative lever 
arms of the plunger gates with respect to the constrictions (see Table[T]). However, no 
correlation is found. G3 is influenced the most by all three plunger gates. G2 and Gi 
are only changed slightly. We therefore conclude that the Coulomb blockade resonances 
and in particular the amplitudes of the current maxima investigated here are mostly 
governed by the wave function in the dot and to a lesser extent by localization sites in 
the leads. 

Furthermore, it can be seen that in general the eflFect of PG3 is opposite to the 
effect of PG2 and PGi. This might indicate that it is indeed not random how the G's 
are changing and that such measurements could be used to obtain further qualitative 
understanding on how the dot wave function is distributed in the dot. The lack of any 
geometric correlation with the evolution of the G's also suggests that the plunger gates 
tune the dot wave function as a whole, rather than several independent puddles. 

4. Conclusion 

We have investigated a three-terminal graphene quantum dot in the multi-level Coulomb 
blockade regime. The dot was thoroughly characterized by both gate-gate sweeps where 
all relative lever arms could be extracted, and temperature dependent measurements. 
When investigating single Coulomb peaks in more detail a shift in peak maxima 
between conductances measured in the different leads were observed. This result can be 
qualitatively reproduced by a rate equation model where different leads couple differently 
strong to different dot states. The shift can be tuned by the plunger gates and by 
investigating the corresponding evolution of the individual conductances we flnd no 
correlation between this evolution and the relative lever arm determined. We therefore 
conclude that on small energy scales the changes in coupling are due to changes in the 
dot wave function, which is rather a single wave function extended over the dot than 
several localized states. This is an important insight in view of the potential to use 
graphene quantum dots for spin qubits. 
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